Dynamics at and near conformal quantum critical points 



S.V. Isakov,' P. Fendley,^'^ A.W.W. Ludwig,"* S. Trebst,^ and M. Troyer' 

' Theoretische Physik, ETH Zurich, 8093 Zurich, Switzerland 
^ Department of Physics, University of Virginia, Charlottesville, VA 22904-4714 USA 
''Microsoft Research, Station Q, University of California, Santa Barbara, CA 93106 
'^Physics Department, University of California, Santa Barbara, California 93106 
(Dated: March 28, 201 1) 

We explore the dynamical behavior at and near a special class of two-dimensional quantum critical points. 
Each is a conformal quantum critical point (CQCP), where in the scaling limit the equal-time correlators are 
those of a two-dimensional conformal field theory. The critical theories include the square-lattice quantum dimer 
model, the quantum Lifshitz theory, and a deformed toric code model. We show that under generic perturbation 
the latter flows toward the ordinary Lorentz-invariant (2-1-1) dimensional Ising critical point, illustrating that 
CQCPs are generically unstable. We exploit a correspondence between the classical and quantum dynamical 
behavior in such systems to perform an extensive numerical study of two lines of CQCPs in a quantum eight- 
vertex model, or equivalently, two coupled deformed toric codes. We find that the dynamical critical exponent z 
remains 2 along the [/(l)-symmetric quantum Lifshitz line, while it continuously varies along the line with only 
Z2 symmetry. This illustrates how two CQCPs can have very different dynamical properties, despite identical 
equal-time ground-state correlators. Our results equally apply to the dynamics of the corresponding purely 
classical models. 

PACS numbers: 64.60.Ht,64.60.Fr,71.1().Pm,71.10.Hf 



I. INTRODUCTION 

The study of critical properties of classical statistical me- 
chanics systems with stochastic relaxational dynamics has a 
venerable history.' More recently, a class of quantum systems 
closely related to classical systems with stochastic dynamics 
has come under intense study. Here, each basis element of the 
Hilbert space of the two-dimensional quantum system corre- 
sponds to a configuration in a two-dimensional classical sys- 
tem. This in itself is not unusual; what is special about this 
class is that the Hamiltonian is chosen so that the ground state 
is written in terms of the Boltzmann weights of the classical 
model. One significant consequence of this relation between 
quantum and classical models is that coiTelators in the ground 
state of such a Hamiltonian are the same as those of the classi- 
cal model. This means that the phase diagrams of the quantum 
and classical models are closely related. 

The quantum dimer model^ is an example of this class of 
quantum systems. On the square lattice, it is a canonical ex- 
ample of the conformal quantum critical points to be discussed 
in this paper, while on the triangular lattice is a canonical ex- 
ample of topological order The Hilbeit space is spanned by 
close-packed hard-core dimer configurations. The quantum 
ground state is the equal-amplitude sum over all such dimer 
configurations. A Hamiltonian having this ground state can 
easily be constructed by using a technique first developed 
in this context by Rokhsar and Kivelson^ (RK). It typically 
is a sum of local projection operators, each annihilating the 



ground state, and will be reviewed in Sec. II B 



In particular, this provides a way to discover and analyze 
new quantum critical points in two spatial dimensions. When- 
ever a ground state in two dimensions has a correlator of lo- 
cal operators algebraically decaying with distance, all local 
Hamiltonians must be gapless."* Thus if the ground state of 
an RK Hamiltonian is described by a critical classical the- 



ory, the quantum theory must be critical as well. Since well- 
understood quantum critical points in two dimensions are few 
and far between, such RK Hamiltonians are quite interesting.!^ 

Indeed some of the best-understood quantum critical points 
in two dimensions were found by analyzing RK Hamiltonians. 
A great deal is known about critical points in rotationally in- 
variant two-dimensional classical models, because these crit- 
ical points are not only scale invariant, but conformally in- 
variant as well. Two-dimensional conformal symmetry has an 
infinite number of generators, and so is very powerful. It can 
be used not only to identify many classes of critical points, 
but to explicitly compute correlation functions in the scaling 
limit, even when the systems are strongly interacting. Because 
the ground state coiTelators are those of the classical theory, 
these results can then be carried over to the quantum case. In 
fact, the ground state itself becomes conformally invariant in 
the scaling limit, as detailed in Ref 7 For this reason, these 
theories were dubbed conformal quantum critical points. 

The static behavior of conformal quantum critical points 
is well understood because of this connection to conformal 
field theory. Their dynamical behavior, however, is another 
story. Only in an exactly solvable casd^^lS dubbed the quan- 
tum Lifshitz theory, has the quantum dynamics has been stud- 
ied in depth. Here the ground-state coiTelators can be writ- 
ten in terms of classical correlators of a free massless scalar 
field. Moreover, the quantum Hamiltonian remains quadratic 
(although not Lorentz-invariant), and so everything in princi- 
ple can be computed.l^l 

One purpose of this paper is to study the dynamics of more 
complex conformal quantum critical points. A key example, 
which we will study in detail, are the two lines of conformal 
quantum critical points in the quantum eight-vertex model of 
Ref 7 The first of these two lines exhibits an additional C/(l) 
symmetry, which reduces the quantum eight-vertex model to 
a quantum six-vertex model whose scaling limit corresponds 
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precisely to the quantum Lifshitz model mentioned above. 
So, while the quantum six-vertex model is not exactly solv- 
able, the field theory describing its scaling limit is. There is, 
however, an RK Hamiltonian whose ground state is that of 
the classical six-vertex model. The second critical line corre- 
sponds to the full quantum eight-vertex model. This line is not 
exactly solvable on the lattice or in the scaling limit, although 
again the ground state of an RK Hamiltonian can be found 
explicitly. The equal-time correlators on the two critical lines 
can be mapped onto each other, because the two correspond- 
ing classical theories are dual to each other.E2l 

In general, different dynamical symmetry classes are possi- 
ble for a given static universality class, depending, e.g., on 
whether the dynamics possesses conservation laws or not.' 
We find here that despite the identical equal-time correlators, 
the dynamical behavior of the two critical lines is indeed quite 
different. In the quantum Lifshitz theory, the dynamical crit- 
ical exponent remains z = 2 all along the line, while along 
the other line, this exponent is found to vary. In the language 
of Ref. [11 the latter possesses Model A dynamics, i.e. the dy- 
namics has no conservation law, and only Ising Z2 symmetry. 
In the former, however, the dynamics respects a U{1) sym- 
metry, and so has a conserved quantity; this corresponds to 
Model B. Both theories possess an exactly marginal operator, 
and the exponents of the equal-time correlators vary continu- 
ously along the critical lines. However, we note that the pres- 
ence of an exactly marginal operator implies the existence of 
a C/(l) symmetry only in two classical dimensions, not in a 
full two-dimensional quantum theory. 

To do this analysis, we exploit the fact that the connection 
between the classical and quantum models goes deeper than 
the equal-time correlators. As described in refs. 1811 1112 ', much 
can be learned about the quantum dynamics by studying the 
classical ones. This means in particular that in some cases 
the dynamical critical exponent z in the quantum model can 
be found by doing classical Monte Carlo simulations. Since 
one point along the z ^ 2 line amounts to doing classical 
Ising dynamics, we make contact with decades of numerical 
studies'^* here. 

Even though an RK Hamiltonian is fine-tuned, when it de- 
scribes a phase with topological order, this physics persists un- 
der (at least small) deformations. This has been demonstrated 
numerically in several examplesj ^ * '^'^ ' and recently been de- 
rived rigorously."^ '"^ The robustness of topological order near 
RK points is not surprising, given that the phase is gapped. 
The argument of course does not apply at conformal quantum 
critical points, and there is no reason to expect that these will 
remain critical under generic perturbations. Even though the 
square-lattice quantum dimer model is critical with seemingly 
no fine tuning, this is a consequence of the highly constrained 
behavior of dimer models. In a more general setting, it has 
been shown that this quantum critical point typically has sev- 
eral relevant perturbations, as well as a host of dangerously 
irrelevant operators.^" ^' 

Thus an interesting question is if a conformal quantum crit- 
ical point is isolated, or if it is part of a phase boundary. In 
renormalization group language, the question is whether a rel- 
evant perturbation causes a flow to another quantum critical 



point, or into a gapped phase. We will present substantial ev- 
idence that the Ising conformal quantum critical point is con- 
tinuously connected to the usual 2H-ld Ising critical point. 

The outline of the remainder of the paper is as follows. In 
Sec. |ll]we review some well known models with RK Hamil- 
tonians, the quantum Lifshitz theory and the quantum dimer 
model. We also display a general connection between clas- 
sical dynamics and quantum dynamics in theories with RK 
Hamiltonians. In Sec. Ill we discuss the dynamics of the 



Ising conformal quantum critical point in two spatial dimen- 
sions, i.e. the quantum critical point whose ground state is 
written in terms of the Boltzmann weights of the critical clas- 
sical 2d Ising model. RK Hamiltonians are strongly fine- 
tuned, and so we show in Sec. [iVjthat generic perturbations of 
the lattice model generate a crossover from the d — 2 dimen- 
sional dynamics with dynamical critical exponent z sa 2.167 
(and 2D static correlation length exponent v — 1), to the criti- 
cal dynamics of the (2H-l)-dimensional classical Ising univer- 
sality class with z = 1 (and ly « 0.632). In Sec. [V] we 
discuss the case of the dynamics of two coupled copies of the 
deformed toric code, or in the equivalent classical case, two 
coupled copies of the critical Ising model. There are two crit- 
ical lines, whose equal-time correlators can be mapped into 
each other's by two-dimensional classical duality. We show 
that these lines correspond to different dynamic universality 
classes. In the case without U{1) symmetry, we will present 
numerical results that indicate that along with the static crit- 
ical exponents, the dynamical critical exponent z also varies 
with the couplings. 



II. CONFORMAL QUANTUM CRITICAL POINTS AND 
RK HAMILTONIANS 

The great progress made in understanding conformal quan- 
tum points came from considering specific ground state wave- 
functions. Once a particular ground state is specified, it is usu- 
ally (but not always) straightforward to find some RK Hamil- 
tonian with this ground state. 

All the models we study have Hilbert spaces whose basis 
elements are labeled by configurations in some classical two- 
dimensional model. Let C label a classical configuration, and 
w{C) be its Boltzmann weight. Then let \C) be the corre- 
sponding basis element of the quantum Hilbert space. Here 
we take the simplest choice for the inner product, the or- 
thonormal one. In a lattice model with discrete degrees of 
freedom. 



{C'\C) = J, 



cc 



(1) 



With continuous degrees of freedom as in a field theory, the 
Kronecker delta is replaced by a Dirac delta function. 

To find a conformal quantum critical point (CQCP) by this 
method, the classical model must be critical and isotropic, and 
the Boltzmann weights w(C) must be real and non-negative 
for all C. Then the (unnormalized) ground state wave function 



\^) = J2VMC)\c) 



(2) 
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is that of a CQCP. The expectation value of any diagonal op- 
erator T) in the ground state is identical to that found in the 
classical theory: 



where V{C) = {C\V\C) is by definition the value of V in 
configuration C. 

At first glance, it appears that for any local dynamics possi- 
ble in a critical classical system, there exists a corresponding 
conformal quantum critical point. This is not true, because in 
a quantum theory one of course weights by the absolute value 
of the amplitude squared. Some classical Euclidean field the- 
ories like the Wess-Zumino-Witten models are only critical 
when the action includes an imaginary term. Thus even if one 
has an RK Hamiltonian whose ground state includes such a 
term, this term disappears from and so the correlators 
will not decay algebraically .- 



a zero-energy ground state. The ground state Q with 



w 



is indeed annihilated by all Q{x). 



Since the quantum Lifshitz Hamiltonian is quadratic in the 
fields, it is exactly solvable, and so much more than just the 
ground-state correlators can be computed.'^ In particular, by 
writing down the three-dimensional Euclidean action for this 
theory 



(f'xdt 



Ik 



(5) 



one sees immediately that the dispersion relation is cx k"^, 
and so the dynamical critical exponent is 2; = 2. As we will 
discuss below, we believe that a generic C/(l)-invariant CQCP 
has z = 2. 



B. Quantum dimers on the square lattice 



A. The quantum Lifshitz theory 

An important example of a CQCP is given by the quantum 
Lifshitz theory. This is an exactly solvable field theory that 
describes, among other things, the scaling limit of the square- 
lattice quantum dimer model. Here the degrees of freedom 
are given by a scalar field (p{x), which is a smooth map of 
all points x on some two-dimensional manifold to a circle. 
The Boltzmann weight of the classical model for a given field 
configuration is simply w — e^'^^'', where S2d is the standard 
action for a free massless scalar field: 

S2d{v) = nj d?x{V^f , (3) 

with K an arbitrary parameter. As discussed in Ref. 7, this 
weight defines a ground-state wave functional, such that all 
diagonal correlators in the quantum model are those of a clas- 
sical massless scalar field in two dimensions. 

The action is invariant under conformal transforma- 
tions of X, so the Boltzmann weights for the quantum Lifshitz 
ground state are invariant as well. It also possesses an addi- 
tional symmetry, under shifts (constant in space) of if. This 
U{1) symmetry is also a symmetry of the quantum Lifshitz 
Hamiltonian 

Hql = IJ d^x [tf + n'iV^f] , (4) 

where the canonical conjugate field n(a;) obeys 
\il{x)^ip{x')] — 16^"^^ [x — x'). For the examples of in- 
terest here, the field is periodic, i.e. ip — tp -\- 2ti. The 
resulting vortex operators, however, are irrelevant at a CQCP.^ 

This Hamiltonian is not a sum of projection operators like 
the RK lattice Hamiltonian is. However, it is very similar: 
it can be written^ as the integral over local terms of the form 

{x)Q{x). Such a Hamiltonian necessarily has non-negative 
eigenvalues, and so any state annihilated by Q{x) for all x is 



The square-lattice quantum dimer model is an example 
of a lattice model with a conformal quantum critical point.l^l 
Each basis element C corresponds to a configuration of dimers 
stretching between neighboring sites on the square lattice. The 
weight wiC) — 1 for all configurations that have exactly one 
dimer touching each site, and is zero for all configurations vi- 
olating the constraint. It has long been known that correlators 
in the classical dimer model on the square lattice are alge- 
braically decaying. Thus by Hastings' theorem,Elthe square- 
lattice quantum dimer model is critical. 

The RK Hamiltonian here was described in the original pa- 
per by Rokhsar and Kivelson.^ It is a sum of local projection 
operators with the ground state annihilated by each projector 
individually. Each term acts non-trivially on a single plaque- 
tte, projecting each of the two configurations with two dimers 
on this plaquette onto their difference. It annihilates all other 
configurations on that plaquette. Each such projector indeed 
annihilates the equal-amplitude sum over all configurations. 

An elegant argument** indicates that the scaling limit of the 
square-lattice quantum dimer model is = 2n of the quan- 
tum Lifshitz model: First, one identifies the classical model as 
having the free scalar field action in ^ with kT^ — 2ti. 
This follows from a Coulomb-gas mapping, done here by 
rewriting the dimer degrees of freedom in terms of a classi- 
cal "height", an integer-value degree of freedom on every site 
of the lattice. It is then natural to assume that the height be- 
comes the continuous field (f in the scaling limit. The value 
K is then identified by comparing the scaling dimensions de- 
termined from the exact computation to those in the classical 
2d scalar field theory. The ground state in the scaling limit 
therefore must be of the form (j2]), because all its diagonal 
correlators are the same as those of S2d- The Hamiltonian 
Q by construction annihilates this. Moreover, note that the 
quantum dimer model has a [/(I) symmetry. The Hamilto- 
nian preserves both the overall number of dimers and their 
"winding number" of dimers. (The winding number is defined 
by nj{—iy, where j labels the links around a cycle and 
rij = 1 for occupied links and rij = for unoccupied.) This 
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FIG. 1 : The finite-size gap of thie quantum dimer model as a function 
of the system size. The gap scales as A oc £,-2.01(2) jjjdicating that 
z = 2.01(2). 



A slight re-writing of the master equation (j6]l gives a gen- 
eralization of the RK Hamiltonian. The rescaled transition 
matrix 



Wc,c' = 



1 



Wc,c' VMC^ 



(9) 



is (real) symmetric due to detailed balance Eq. ([8]l. Letting 
PC (t) = PC (t) / ^/w{C), the rewritten master equation is 



dT 



PC 



iT) = J2 Wc^c'PC'ir) 



(10) 



This is a Schrodinger equation in imaginary time r with 

Hamiltonian Hc.c — l)VFc,c' for the time-dependent 
quantum state 



(11) 



requires that only derivatives of ip can appear in the Hamilto- 
nian, which indeed is a property of Q. 

Finding the dynamical scaling exponent z provides another 
check that the scaling limit of the square-lattice quantum 
dimer model is given by the quantum Lifshitz theory. The 
quantum dimer model is not solvable, even though its ground 
state is known exactly. We thus need to find z numerically; 
we do this by studying the finite-size scaling of the gap. Our 
method is described in Appendix[A| We find that z — 2.01(2), 
as shown in Fig.[T] This is in excellent agreement with the ex- 
act value z = 2 for the quantum Lifshitz theory, as seen from 
(jsjl, and previous numerical results.^- As such it also provides 
a nice consistency check on our numerical method. 



C. RK Hamiltonians from classical stochastic dynamics 



For classical systems with (positive) Boltzmann weights 
w{C), every relaxational stochastic dynamics is known" to 
give rise to a "RK Hamiltonian": time-dependent probabili- 
ties pc{t) > 0, satisfying a "master equation". The master 
equation is written in terms of a "transition matrix" Wc.c, 
where for different configurations C ^ C, the matrix element 
Wcfi' > denotes the transition rate from C to C. The di- 
agonal element defined by Wc,c '■= — Tlc'^c ^C'.c denotes 
the transition rate out of state C. The master equation is then 



c 



Wc,c' Pc'{t) 



(6) 



The probabilities relax at long times to the classical equilib- 
rium distribution 



Pc{t) w{C) 

when detailed balance 

wiC) Wc'fi = Wcfi' w{C') 

is satisfied. 



(7) 



(8) 



Note that Wc,c' and Wcc are related by a similarity transfor- 
mation and so have the same spectrum. The Hamiltonian thus 
has the same spectrum as that of the relaxational dynamics, 
Eq. (j6|l. The state ( [TT] i relaxes, due to Eq. (j7]i, at long imagi- 
nary times T to the ground state Eq. (j2]). Thus this construction 
indeed provides a generalization of the RK-type Hamiltonian 
and ground state. 

In the present formulation, C denotes configurations of clas- 
sical variables on a lattice. When the system with classi- 
cal Boltzmann weights w(C) possesses an equilibrium critical 
point, a universal continuum field theory description emerges 
upon coarse graining for both the static as well as for the 
dynamic correlations. The resulting dynamic universality 
classes are classified within the well known framework of 
Hohenberg and Halperin. ' For example, the Ising conformal 



quantum critical point discussed below in section 111 is de- 



scribed by the Model A dynamical universality class of Ho- 
henberg and HalperirP (describing dynamics lacking any con- 
servation law). The quantum Lifshitz universality class, how- 
ever, possesses a U{1) symmetry, so that this belongs to the 
Model B universality class. 

The classical dynamics are conventionally described by 
a time-dependent Landau-Ginzburg equation with stochas- 
tic Langevin-type noise, or equivalently, as a Fokker-Planck 
equation for the time-evolution of the probability distribu- 
tion for the coarse-grained classical degrees of freedom. It is 
well knowrP^ that if one performs the corresponding similar- 
ity transformation to make the linear operator appearing in the 
Fokker-Planck equation Hermitian, the Fokker-Planck equa- 
tion turns into the coarse-grained analog of the Schrodinger 
equation, Eq. (lOi, in imaginary time. Specifically, for a 



classical statistical mechanics system in d spatial dimensions 
described, e.g., by a static (real) Landau-Ginzburg action 
S{4)a{x)} for coarse-grained classical degrees of freedom 
{4>a{x)}, the resulting quantum Hamiltonian is again of the 
form 



H = 



d^'x Ql{x)Qa{x) 



(12) 
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generalizing Eq. d?]), where 



Qa{x) 



6 ^1 SS 

d(l)a{x) 2S(f)aix) 

6 1 SS 



1 

71 
1 



(13) 



This implies very generally that the ground state wave func- 
tion ^{{(j)aix)}) = {l/VZ)exp{-lS{(l)a{x)}} is given by 
the classical Boltzmann weight (Z denotes the classical parti- 
tion function). Thus the (in general non-Hermitian) linear op- 
erator appearing in the Fokker-Planck equation becomes, after 
similarity transform, the negative of the (Hermitian) quantum 
Hamiltonian, {—1)H. All its eigenfunctions are positive. 

The Hamiltonians of the form of Eq. (12 1 are clearly very 
special. One may ask if there exists a symmetry principle 
which constrains Hamiltonians to be fine-tuned to this very 
particular form. Indeed, there exists such a symmetry. It is 
again well known^"* that classical relaxational stochastic dy- 
namics can be viewed as being invariant under certain su- 
persymmetry transformation.^^ The property that equal-time 
correlation functions converge at large times to static equilib- 
rium correlation function can then be seen as a direct conse- 
quence of the underlying supersymmetry. The perturbation of 
the Ising conformal quantum critical point discussed in Sec. 
IV can thus be viewed as breaking the supersymmetry. 

This procedure makes obvious the relation between classi- 
cal dynamics and RK Hamiltonians. In fact, the generic way 
of constructing RK Hamiltonians in field theory discussed in 
Ref is formally identical to that in Eqs. ( [T2] l and ( [T3| . For 
example, the quantum Lifshitz Hamiltonian (pli arises with 
one field ip and setting S — S2d from Eq. ([sjl.'^l^ 



III. QUANTUM AND CLASSICAL DYNAMICS IN THE 
ISING MODEL 



The Hilbert spaces of the lattice models we consider for the 
remainder of the paper are simply those of the Ising model 
on the square lattice. Namely, at every site we have a spin 
(i.e. a two-state quantum system), and all the Hamiltonians 
we consider are invariant under flipping all the spins. We write 
the Hamiltonians in terms of the Pauli matrices acting on 
the spin at site a. 

For both pedagogical and physical reasons, it is often con- 
venient to describe the Hilbert space in terms of domain walls 
instead of spins. There is a two-to-one mapping from spin 
configurations to domain walls obtained simply by drawing 
a lines on the links of the dual lattice separating regions of 
different spins. A typical configuration is illustrated in Fig. [2] 
Each domain wall is a two-state system, so we define the Pauli 
matrices ct^ acting on the link k analogously to the spins. 
However, the domain walls cannot end or branch, so there 
must be an even number of domain walls at each site. In an 
equation, we constrain 



n 



= 1 



(14) 
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FIG. 2: Sketch of a typical domain wall configuration separating re- 
gions of opposite spins. 



for all sites a on the dual lattice. For the obvious reason, these 
domain walls are often referred to as "loops". 



A. The toric-code ground state 

Let us first consider an important non-critical case, closely 
related to the mother of all lattice models for topological quan- 
tum computation, the toric code.^^ Here the ground state ^ 
is a sum over all Ising configurations with the same ampli- 
tude w{C) — 1. Thus the associated classical model is sim- 
ply the Ising model at infinite temperature. Equivalently, the 
ground state can be viewed as a sum over all loop configura- 
tions with equal amplitude. Correlators of local objects in the 
ground state (e.g. whether a link is occupied by a loop or not) 
are obviously finite-range. 

The gapped RK Hamiltonian having this ground state is 
simply the sum over all possible single-spin flips: 



(15) 



In terms of domain walls/loops, this acts on the four links on 
the dual lattice surrounding each site j, removing a wall when 
there is one and adding a wall when there isn't: 



n 



(16) 



k touching 



where the product is over the four links on the plaquette on 
the dual lattice surrounding the site a. 

This equal-amplitude sum over all loops/domain walls is 
the same ground state as in the toric code. The difference be- 
tween the models is that in the toric code, the constraint ^\A\ 
is not imposed. Rather, a diagonal term is added to the Hamil- 
tonian to give an energy for every sites with an odd number of 
domain walls touching it. The energy penalty then precludes 
any configuration with an odd number of domain walls at a 
site from appearing in the ground state. 

Since each term in the Hamiltonian ( [Ts] ) commutes with 
each of the others, this model is trivially solvable, and the 
sum over all loops is indeed the ground state. On the disk, this 
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ground state is unique. On the torus, there are four ground 
states, corresponding to the choices of having zero or one 
loops go around each cycle. It is easy to see that the model 
is gapped; excited states correspond to having plaquettes with 
the eigenvalue —1 of r^. This ground state has topological 
order, and so has been the subject of intense study in recent 
years. 



B. The Ising CQCP 

Thinking about the toric code in terms of Ising spins makes 
it clear how to find a conformal quantum critical point sepa- 
rating the phase with topological order from an ordered phase. 
We simply need to change the weight in the ground state 
to correspond to that in the classical two-dimensional Ising 
model at non-infinite temperature. In loop language, this cor- 
responds to adding a weight per unit length 2ry, so that 

w{L) = jf"" 

where nc is the number of links covered by loop. In the quan- 
tum theory, we can write 



f N \ 

Y,ol\\C)^(2nc^N)\C) 

\k=l ) 



for N links. The ground state is then simply 
|vl/)=^^"(^)|£). 



(17) 



An RK Hamiltonian with the ground state ( 17 1 for general 
77 can be found by including a potential that weights links. We 
will show below that the Hamiltonian found in Ref. [T] gives 
two copies; a simpler-to- write form iJ^H 



(18) 



This Hamiltonian can be written as a sum over projectors, 
each projector breaking into two-by-two blocks. Note that 
there is no manifest C/(l) symmetry like there is in the quan- 
tum dimer model. 

The 2d Ising critical point is at = r\c, where 




r;, = (l + ^/2)-l/2 



(19) 



Any local Hamiltonian such as ( [TS] ) having (17 1 with -q = 
77c as a ground state must be quantum critical by Hastings' 
theorem. We call the Hamiltonian ( [T8] l with 77 — j]c the Ising 
CQCP. 

This quantum Hamiltonian imposes a dynamics on the 
spins identical to the usual classical dynamics with local up- 
dates. Thus the dynamical exponent z is the same in two 
cases. There have been numerical determinations of z for 
decades; the current value from Ref. |29]is z = 2.1667(5). 
This is slightly larger than the z = 2 value for the quantum 
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FIG. 3: The finite-size gap of Hamiltonian ijTS} at the Ising CQCP r/c 
as a function of the system size. The gap scales as A oc 7^-2.170(3) 
indicating that z = 2.170(3). 



dimer model and the quantum Lifshitz theory; below in Sec. 
|V]we explain how to adapt an argument of Ref. 7 to suggest 
that z > 2 for CQCPs. For the Ising case here, a rigorous 
inequality requires that z > 7/4.E2I 

We have rechecked this number using our numerical meth- 
ods described in Appendix [A] and find the same result, see 
Fig. [3] We note that one needs to simulate the classical Ising 
model with the Monte Carlo transition matrix that is propor- 
tional to the quantum Hamiltonian. In this case, the transition 
matrix has the following form 



Wc,c' = V 



lw{C') 



where 'w{C) and w{C') are the weights of the configurations 
before and after update (single spin flip). This transition ma- 
trix leads to the same dynamical exponent as the Metropolis 
transition matrix. 



IV. FLOW FROM 2D ISING TO 3D ISING 

Another interesting way to deform the toric code is to in- 
clude the more conventional Ising nearest-neighbor interac- 
tion. In this section, we study the phase diagram of the toric 
code deformed by this interaction and the special RK-type in- 



teraction in ( 18 1. We will show that including this term causes 
a flow from the Ising CQCP with 2d Ising critical exponents 
to the usual z = 1 Ising critical transition in 2 + 1 dimensions. 
Adding the Ising nearest-neighbor interaction deforms 



Eq. (18lto 



H 



-y 

774 



(20) 

where the sites /3 are the nearest neighbors of a. We include a 
factor of 1/77^ compared to Eq. 18 for convenience. In terms 
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FIG. 4: Phase diagram in the h-l/rj plane. Error bars are smaller 
than the symbol size. The line guides the eye. The whole transition 
line has z = 1 except for one point at h = and 1/rj — {1 + 
\/2)^^^. For the two points indicated as A and B, we extract the 
critical exponent from data collapses of the Binder cumulants in 
Figs. |5] and |7] respectively. 



of the link degrees of freedom, the extra term is simply a loop 
tension —h J2p <^p- The other diagonal term in H can also be 
thought of as a loop tension, but fine-tuned to make the ground 
state have the weighting of the 2d classical Ising model. 

The Hamiltonian ( |20| reduces to that of the 2D transverse 
field Ising model on the square lattice for i] = 1. There is a 
phase transition at w 0.32847,"" which is in the 3D Ising 
universality class with the dynamical critical exponent z = 1 
and the correlation length exponent sa 0.632. From the dis- 
cussion in the previous section, we know that there is a confor- 
mal critical point at /i = and 77 — rjc. This phase transition 
is in the 2D Ising universality class with the dynamical criti- 
cal exponent z « 2.167 and the correlation length exponent 
1/ = 1. 



A. Phase diagram 

The interesting question is if these two critical points are 
connected by a line of phase transitions, or describe discon- 
nected regions in parameter space. We find convincing nu- 
merical evidence that the former is true. The full phase dia- 
gram in the h — 1/r/ plane is shown in Fig. |4] The CQCP at 
h = and i] = rjc is unstable under perturbations by both 
h and 77. Along a particular line in this plane, there is a flow 
from the Ising CQCP to the 2H-ld Ising fixed point. The whole 
phase boundary in Fig. |4] has 3D Ising exponents except for 
one point. 

This phase diagram was mapped out by using a variant of 
the continuous-time quantum Monte Carlo algorithm. We 
measure the Binder cumulant in Monte Carlo simulations 
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FIG. 5: Data collapse of the Binder cumulant ([2TJ for point A at /ic = 
0.076, 1/rj = 1.3 using the 3D Ising exponent u = 0.632 (z = 1). 
Data is obtained from continuous-time Monte Carlo simulations at 
inverse temperature /3 = 2L. 



scales in the vicinity of a continuous phase transition as 

U{L,K,f3)^F{L'/''{K~K,),/3/L^), 

where F is the scaling function, L is the linear system size, z 
is the dynamical critical exponent, ly is the correlation length 
exponent, K — Kc is the distance to the critical point in some 
coupling constant K, and f3 is the inverse temperature. It fol- 
lows from the above equation that the curves for different sys- 
tem sizes should collapse onto the universal curve F for ap- 
propriate values of and Kc when is fixed. We locate 
the critical points shown in Fig. |4]by collapsing the Binder 
cumulant data. 

In Fig. |5] we show an example of such data collapse for 
he — 0.076, l/f] — 1.3. It is very hard to obtain the phase 
boundary close to the Ising CQCP using the continuous-time 
algorithm because the gap becomes very small and one needs 
to simulate very low temperatures to obtain any meaningful 
results. We believe that the phase transition is in the 3D Ising 
universality class with z — 1 and v — 0.632 at any point on 
the phase boundary for finite values of h. 



B. Instability of the 2D Ising point to Trotter errors 

In this subsection, we show that the Ising CQCP is un- 
stable even to the Trotter discretization eiTors. The imagi- 
nary time direction can be discretized in iV = /3/e slices 
and using the Suzuki-Trotter decomposition, one maps the 
2D quantum model onto the 3D (Ising-type) classical model 
with inter-layer couplings given by the diagonal couplings in 



where m is the magnetization density. The Binder cumulant 



Eq. (20 1 and the intra-layer ferromagnetic coupling K'^ = 
— 1/2 In tanh eF, where F = 1/77"* is the strength of the quan- 
tum transverse field. K'^ diverges in the limit e — > 0. 

This approach is often used to obtain the critical expo- 
nents when simulations using the continuous-time algorithm 
are cumbersome. The exponents obtained by discrete-time 
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FIG. 6: Data collapse of the Binder cumulant \2l\ for h — and (l/7?)c = 1.316, obtained from discrete-time Monte Carlo simulations with 
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and continuous-time methods should have the same values for 
stable fixed points. However, in our case, the 2D Ising point 
is unstable and we obtain the 3D Ising exponents in discrete 
time indicating that there is a flow to the 3D Ising fixed point, 
see Fig. |6] Thus we may have the case when quantum-to- 
classical mapping fails. In principle, one should recover 2D 
exponents in the limit e — > 0. We are unable to do this because 
the third direction becomes very large for small e and Monte 
Carlo simulations become quite impractical. 



C. PIGS simulations 

The continuous-time quantum Monte Carlo simulations 
around the Ising CQCP become very slow because we use 
local updates in the spatial direction. To prove that the tran- 
sition at /i = and rj = rjc is indeed in the 2D Ising uni- 



versality class, we perform quantum Monte Carlo simulations 
by using the path integral ground state (PIGS) Monte Carlo 
algorithm.^^^ In the PIGS algorithm, one uses a variational 
wave function and tries to project the ground state wave func- 
tion. In our case, the ground state is known exactly (strictly 
speaking, the wave function given in Eq. (j2]i is not the ex- 
act ground state in discrete time but becomes one in the limit 
e — )■ 0). Thus simulations can be performed at any temper- 
ature if one uses the continuous-time algorithm. We choose 
1/T = 0.4. This temperature is low enough to have sub- 
stantial quantum dynamics and to avoid trivial classical sim- 
ulations of the Ising model that one effectively has at high 
temperatures. Figure. [7] shows the data collapse of the Binder 
cumulant ( pT| with the 2D Ising correlation length exponent. 
The 2D Ising exponents can also be recovered using the PIGS 
Monte Carlo algorithm in discrete time at small enough values 
of e, not shown. 
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Carlo simulations at P — 0.4. 



CONFORMAL QUANTUM CRITICAL LINES WITH 
CONSTANT AND VARYING z 



The two-dimensional classical Ising model at its critical 
point provides a simple example of a conformal field theory. 
All the scaling dimensions can be determined, and one finds 
that there are no marginal perturbations possible. However, by 
coupling two critical Ising models together, one finds a theory 
with an exactly marginal operator, leading to a line of renor- 
malization group (RG) fixed points. If the Ising spins are on 
the same lattice, this is called the Ashkin-Teller model, while 
if they are on interpenetrating square lattices, this is called the 
eight-vertex model.tSl In equivalent fermionic language, this 
amounts to coupling two Majorana fermions together with an 
exactly marginal four-fermion term. 

This richer behavior persists in quantum theories with RK 
Hamiltonians. We review below how this coupling yields two 
conformal quantum critical lines. ^ For any given point on one 
line one can map its correlators - given in terms of the under- 
lying classical theory - to those of a point on the other line. 
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FIG. 8: The four configurations in the coupled deformed toric codes, 
and their weights in the ground state. 



Thus the equal-time correlators in the quantum theory have 
the same behavior at these two related points. This, however, 
is no guarantee that the quantum theories have the same dy- 
namical behavior Indeed we will show in the following that 
the dynamical behavior turns out to be quite different. 

On one of these lines, the quantum six-vertex model, the 
model has a height description and so has a [/(I) symmetry. 
Thus one expects the quantum Lifshitz theory to describe the 
scaling limit, and indeed we see numerically that z remains 
2 all along this line. The second conformal quantum critical 
line includes the case where the two Ising CQCPs are decou- 
pled. Thus there is no U{\) symmetry, and as we saw above, 
z is not 2. However, these two critical lines intersect, and so 
necessarily z cannot remain at the Ising value along this line. 
We present numerical work that indicates that z indeed varies 
continuously along this line. 

In fact, an argument adapted from Ref. |7] suggests that 
z > 2 for a CQCP. By relating the stress tensor of the 
two-dimensional conformal field theory describing the ground 
state of the CQCP to deformations of the quantum Hamilto- 
nian, it is argued that terms in the latter can only depend on 
space via the Laplacian, e.g. terms like {V'^LpY in the quan- 
tum Lifshitz theory. In general there is no scalar field, but the 
argument still suggests that the usual Lorentz-invariant terms 
in the effective action yielding z ~ 1 are absent. Moreover, 
fields in a unitary conformal field theory necessarily have non- 
negative dimension, so acting with the dimension-2 Laplacian 
gives a field of dimension at least 2. This suggests that z > 2 
for CQCPs where the underlying conformal field theory is uni- 
tary, as for the theories considered here. It would be interest- 
ing to make this argument more precise; we indeed find that 
this bound is satisfied for all the CQCPs we study. 



A. The model 



The phase diagram of the coupled model is well under- 
stood, because the corresponding classical model is equiva- 
lent to the eight- vertex model This can be seen by first 
reverting to the original Ising spins on the doubled lattice, and 
then drawing the domain walls on the dual lattice to the dou- 
bled lattice. The four configurations in Fig. |8]then become 
the eight vertices of the zero-field eight-vertex model. The 
counting is as follows. There are of course 16 possible con- 
figurations of four Ising spins around a plaquette. One can flip 
all the spins on the lattice of solid circles without changing the 
configurations in Fig. [8] One can do the same for the lattice 
of open circles. Flipping the spins on both lattices does not 
change the eight-vertex configuration, but flipping the spins 
on say the solid lattice does. Thus the 16 different Ising spin 
configurations correspond to the eight different configurations 
in the eight- vertex model and the four different configurations 
in Fig. |8] In the standard eight-vertex model language, the 
four configurations in Figure |8]have weights c, a, h and d, re- 
spectively. To make the connection with coupled Ising do- 
main walls more apparent, we set c = 1 /yy and d = rjx and 
a = 6 = 1. 

The RK Hamiltonian with this ground state is therefore the 
same as that of Ref. H It is comprised of a sum over 2 by 2 
blocks just like the Ising case. The off-diagonal terms flip all 
the solid lines around a plaquette on the dual lattice, or flip the 
dotted lines around a plaquette on the original lattice. The 2 
by 2 block thus acts on four faces surrounding each point on 
the doubled lattice. Let be the number of these four faces 
which are empty, and rid be the number of crossings in these 
four Likewise, let Tic and Ud be the number of empty faces 
and crossings in the flipped configuration. Then the diagonal 
term for each site is given by c^"<^d^"'', i.e. what enters is the 
number of crossings and empty faces in the flipped configura- 
tion. Each 2 by 2 block is then 



Q^ 



orf2n^ 



(22) 



where c ~ l/rj and d = ijx as before. The reason for the 
prefactor as compared to Ref. |7] is to ensure that the terms 
remain finite in the c — > and d limits. The fact that 
nc + nd + ric + nd = 4 can be used to rewrite the off-diagonal 
terms if desired. 



As before, we study the Ising models in their domain 
wall/loop formulation, so that the degrees of freedom live on 
the links of the lattice. We place these models on two inter- 
penetrating square lattices (i.e. on both a square lattice and its 
dual). Then the quantum model has a four-state system on 
each face of the doubled lattice, illustrated in Figure |8] The 
sites on the original lattice are labeled by solid circles, while 
the dual lattice sites are labeled by the open circles, so the 
figures applies to half the faces; the others are given by a rota- 
tion. The constraint ([T4| is required on all sites on the doubled 
lattice, assuring that neither type of loop branches or ends. 



B. Constant z = 2 along the (7(1) symmetric Hne 

One critical line occurs when there is a C/(l) symmetry, on 
the quantum six- vertex line. This U{1) symmetry arises when 
X = 0, i.e. the dotted and solid lines are never permitted to 
cross in the ground state, so that there is a three-state system 
on each face. The symmetry becomes apparent by rewriting 
the degrees of freedom in terms of an integer-valued "height" 
h at each site of doubled lattice. Heights on adjacent sites dif- 
fer by one. The solid lines in Fig. |8]represent domain walls 
between heights on the original lattice, while the dashed lines 
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are domain walls between heights on the dual lattice. On a 
disc or sphere, this definition is unique up to an overall shift 
h ^ h + noia sign flip h — > —h. (When x ^ 0, one can con- 
sistently assign only Z2-valued heights, i.e. Ising variables, 
consistently.) 

Thus this height has the same properties as the field (p in 
the quantum Lifshitz model discussed in Sec. II A and it is 



natural to identify h with tp in the continuum limit. However, 
as with the XY model, the classical model is not critical for all 
values of 77; a Kosterlitz-Thouless transition occurs at c = 2. 
Thus for X — 0, the model is critical for any ?/ > l/\/2. 
This critical line should be described by the quantum Lifshitz 
model. The exponents in equal-time correlators will depend 
on rj, but we should have z = 2 all along this line. We have 
checked this explicitly for one point at x = and 77 = 1 /\/2- 
A subtlety in taking the a; = limit is that crossings are not 
forbidden with this Hamiltonian. Rather, they become fixed 
defects. However, except in peculiar special cases all such de- 
fects have a non-zero gap, and so can be ignored in the scal- 
ing limit. Another thing to note is that there are extra ground 
states on the torus in the a; = limit. These are the analog 
of tilted states in the dimer case, where the height is not pe- 
riodic around a cycle of the torus. In the language of dashed 
and solid loops, these result from configurations where the 
loops around a cycle alternate between dashed and solid. Two 
neighboring loops of the same type around a cycle can anni- 
hilate, but loops of different type cannot: they lie on different 
sets of links. Thus when crossing is forbidden, there is no way 
for non-contractible loops of alternating type to annihilate. 



C. Continuous variation of 2 along tlie Z2 -symmetric line 

The second critical line^ in the quantum eight- vertex model, 
which does not have a [/(I) symmetry, is parametrized by 



± 2, or 



(23) 



We will refer to this line as the "Z2 critical line". For Xc = 
this line meets the U{1) critical line discussed in the previ- 
ous section, and we thus expect a dynamical critical exponent 
z ~ 2 at this point. However, when 77^ = \/2 — 1 on the 
Z2 line, we recover the case of two decoupled Ising models 
corresponding to — 1. For this second point we thus ex- 
pect a dynamical critical exponent of z w 2.17 as previously 





FIG. 9: Two configurations differ by a single plaquette flip. Thin 
lines denote empty faces and thick lines denote nonempty faces. 
Here ric = 2, = 0, = 0, and hd — 2, see text. 



obtained in extensive numerical simulations of the dynamics 
of the 2d Ising model,'^ see also Fig. [3] As a result, one is 
led to expect that the dynamical critical exponent is varying 
(continuously) along the Z2 critical line. This might not be 
surprising since it is known that all static critical exponents 
(of the equal-time correlators) are continuously varying along 
the Z2 and U{1) critical lines. Thus it might be natural to ex- 
pect that in the absence of additional symmetries all critical 
exponents, including z, are continuously varying along such a 
critical line. 

To calculate the variation of the dynamical critical expo- 
nent along the Z2 critical line, we again perform classical 
Monte Carlo simulations. The off-diagonal part of the Hamil- 
tonian of the two coupled deformed toric codes is more com- 
plicated compared with the single deformed toric code. To 
work out the transition matrix, we note that the weights be- 
fore and after update are w{C) = (^xr^Y^'^rf''^''^'^^'^ and 

respectively. An example of 



2\2hd 2nc_+2nd 



w{C') = {xri'^y"'''q' 
such an update is shown in Fig. l9] 
matrix element Qij = —(1/77)"'^^'' 
obtain the transition matrix 



Given the off-diagonal 

(a;77)"''+"'' , it is easy to 



Wc,c' = w{C'). 
This is indeed a legitimate transition matrix since < w < 1. 



In Fig. 10 we show our numerical estimates of the dynam- 
ical critical exponent z as a function of the coupling x. The 
point a; = is the special point where the two critical lines 
meet. Here z = 2 (within error bars), as it remains all along 
the critical line with U{1) symmetry (x ~ and 77 > l\/2)- 
Moving along the Z2 critical line with increasing coupling x, 
we see that z indeed varies continuously. The dynamical ex- 
ponent has a non-monotonic behavior with an increase up to a 
maximum of 2.196 close to the Ising point (x — 1) and a de- 
crease down towards z = 2 as x goes to infinity. We note that 
there might be small systematic errors, which are discussed in 
Appendix [B| in detail, because the system sizes are not large 
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FIG. 10: The dynamical exponent z as a function of x along the crit- 
ical line. The dashed line indicates the coupling x = 1 of the case of 
two decoupled classical Ising models. Error bars shown correspond 
to the statistical errors of the numerical fit described in Appendix [A| 
For a discussion of systematic errors see Appendix [B| 
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enough for the majority of points in Fig. 10 

The contrasting observation that the dynamical critical ex- 
ponent does not change along the U{1) critical line - despite 
continuously varying equal-time correlators - suggests that 
this additional U{1) symmetry protects the dynamical critical 
exponent z — 2. This is in line with the general classification 
of Hohenberg and Halperin ^ within which the Z2 symmet- 
ric line is in the Model A dynamical symmetry class (absence 
of any conservation law), and the U{1) symmetric line in the 
Model B dynamical symmetry class (presence of a conserva- 
tion law). 



VI. DISCUSSION 

We have studied the dynamics of a number of conformal 
quantum critical points. We explored the Ising CQCP in 
depth, showing that its dynamics are equivalent to classical 
dynamics. We also showed that it is quite unstable, and that 
by perturbing with the usual nearest-neighbor interaction it 
can be continuously connected to the conventional Lorentz- 
invariant Ising transition in 2+\ dimensions. By studying two 
coupled deformed toric codes, we illustrated how different 
dynamics can result in different values of z, even when the 
equal-time correlators are the same. 

There are several interesting directions for future research. 
At the beginning of Sec. [V]we discussed two quantum-critical 
lines, one possessing U{1) symmetry and the other possessing 
only Z2 symmetry. These two lines intersect at one point at 
which the dynamical critical exponent is z = 2. This point, 
not surprisingly, has enhanced symmetry. For example, the 
equal-time correlations are known to possess an SU{2) sym- 
metry. Moreover, the physics should be in the same univer- 
sality class as the loop models studied in Ref. [34] The cou- 
pling constant moving one along the Z2 line away from this 
intersection point with SU (2) symmetry must be an exactly 
marginal perturbation which generates a line of fixed points 
for both statics and dynamics. In particular, by perform- 
ing perturbation theory in this coupling constant around the 
SU{2) point, one expects to be able to obtain analytic expres- 
sions for the deviation of the dynamic critical exponent from 
z = 2. 

It would be interesting to explore whether the quantum crit- 
ical line recently postulated'^ for the toric code model in a 
multicomponent magnetic field bears some relation to our re- 
sults. It was found'^ that along this line the product zv of 
dynamical critical exponent and correlation exponent appears 
to vary from approximately 0.69 to 1. The most likely in- 
terpretation of those results might be in terms of a crossover 
between a conformal quantum critical point with z — 2 and 
correlation exponent v — 1 /2 - corresponding to the x 00 
limit of the Z2 critical line in the quantum eight-vertex model 
studied in this manuscript - and a Lorentz-invariant (2h-1) di- 
mensional multicritical point. Such a scenario would be akin 
to the crossover discussed in detail in Sec. |IV] 

One may further ask if it is possible to define a quantum 
dynamics of RK type based on classical stochastic dynamics 
as discussed in this paper, for (2h-1) dimensional systems sup- 



porting non-Abelian statistics. A case in point is the Levin- 
Wen model, ■'^ which can be viewed-'^ •'^ as a 'lattice regular- 
ization' of (2H-l)-dimensional doubled SU(2) Chem-Simons 
theory at level fc = 3, possessing anyon excitations with non- 
Abelian exchange statistics. Indeed, the ground state j^*) of 
this model in the simplest non-abelian (Fibonacci) case can 
be described in a geometric form similar to Eq. Such a 
description is found by using the results of Ref 38 to rewrite 
the Fibonacci Levin- Wen model in terms of the quantum net 
model of Refs. '39'40' Here, the configurations C describe 
configurations of so-called "nets", which provide an orthonor- 
mal basis as in Eq. ([TJ. A key difference with Eq. (|2]l is that 
w{C) is replaced by a real wavefunction of the configura- 
tions C which takes on both positive and negative values. (The 
lack of positivity of the wavefunction ("sign problem") may 
be a general feature of systems supporting excitations with 
non-Abelian statistics.) This property prevents one from con- 
structing a quantum dynamics of RK type for the Fibonacci 
Levin- Wen model by using a stochastic relaxational dynamics 
of a suitable classical statistical model, as we did in Sec. |IIC| 

Nevertheless, there is a classical 2D statistical mechanics 
model which arises naturally from the ground state \^) of the 
Fibonacci Levin- Wen model, and gives rise to a CQCR^^^ 
By tuning the weight per unit length of net just as we did 
in this manuscript, the (2H-0)-dimensional classical partition 
function obtained from the wave function j'l'), Z — (^'|^) is 
critical. Z is now a sum of non-negative Boltzmann weights 
because only the square of the wavefunction appears due to 
the orthonormality of the basis. By Hastings' theorem, the 
quantum model must therefore be critical. This classical criti- 
cal point is in the universality class of the (2 + 0)-dimensional 
conformal minimal model with central charge c = 14/15 
(or in the language of Ref. I39I40I the Q-state Potts model 
with Q = (5 + \/5)/2).'^ Such a conformal field theory is 
knowri^to be described by a classical Boltzmann weight aris- 
ing from Landau-Ginzburg theory, precisely of the type de- 
scribed above Eq. ( [T2j l. Were it not for the "sign problem" 
mentioned above, one could proceed to study the dynamic 
critical exponent z for this conformal quantum critical point 
along the lines of the present manuscript. 
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FIG. 11: Four times the gap of Hamiltonian \IS\ and the gap of 
Hamiltonian \22\ {x — 1) at the critical point rjc as a function of 
the system size. The gaps are multiplied by L^'^^ to make the dif- 
ference between the two curves clearly visible. The gaps of the two 
decoupled Ising models show very strong system size dependence. 
Fitting the curve from L = 24 to L = 64 gives z — 2.192(3) (indi- 
cated by the dashed line), whereas fitting the curve from L = 64 to 
1/ = 96 gives z — 2.170(4) (indicated by the solid line) as for the 
single Ising model. 



Appendix A: Numerical calculation of the dynamical critical 
exponent 

In this Appendix, we describe our numerical method. To 
extract the dynamical critical exponent z, we use the finite 
size scaling of the gap A oc at the critical point, where 
L is the linear system size. The gap can be obtained from 
the autocorrelation function of classical Monte Carlo simu- 
lations by ensuring that the off-diagonal part of the Monte 
Carlo transition matrix is proporti onal to the off-diagonal part 
of a Hamiltonian** (see also Sec. II C I. Strictly speaking, the 
Hamiltonian is symmetric but the transition matrix Wc,C' is 
not. The following matrix is symmetric and has the same 
eigenvalues as Wc,C'-^ 



e.c 



w{Cf/^Wc,c'w{Cr"\ 



where w{C) and w{C') are the statistical weights of the con- 
figurations before and after update. If the transition matrix W 
is proportional to the Hamiltonian: 

W = l-cH, 



where c is the coefficient of proportionality, then the gap of 
the Hamiltonian is related to the autocorrelation time ta of 
some observable A as 



1 



A 



where ta is measured in Monte Carlo time units. The autocor- 
relation time Ta can be obtained by fitting the autocorrelation 
function to an exponential function (at long enough Monte 
Carlo time to make sure that the contribution from the other 
modes is negligible). 

We perform Monte Carlo simulations as follows. We sim- 
ulate a classical model corresponding to the quantum Hamil- 
tonian for different system sizes L. Typically the largest sys- 
tem size is L = 64. We measure the autocorrelation time of 
the magnetization for Ising-Iike models and the autocorrela- 
tion time of a product of two dimers separated by a distance 
of L/2 for the dimer model. The dynamical exponent is cal- 
culated as described in the previous paragraph. Typically we 
perform 25 independent runs for every point and the error bars 
are obtained by using the jackknife method. 



Appendix B: Size dependence of the gap in the deformed toric 
code model 



We note that at the decoupling point x = 1 of the deformed 
toric code model, describing two decoupled Ising models, the 
gap shows very strong system size dependence reaching the 
asymptotic exponential form only at large system sizes. In 
Fig.[TT[ we show the gaps of this model and the ordinary Ising 
model as functions of the system size. The dynamical critical 
exponent of the two decoupled Ising models is z = 2.192(3) 
if we fit the curve from L = 24 to L = 64 and z = 2.170(4) if 
we fit the curve from L = 64 to L = 96. The latter value is in 
agreement with that for the ordinary Ising model.l22l The ori- 
gin of this strong system size dependence might be tracked 
back to the additional prefactor c"'^+"<=c?"''+"'' in Hamilto- 



nian (22 1 as compared to the ordinary Ising Hamiltonian ( 18 1. 
This prefactor might give rise to slightly modified dynamics 
for small system sizes. Typically the largest system size in 
our Monte Carlo simulations is L = 64 for small values of x 
and i = 48 for large values of x. Assuming that the gaps for 



all points in Fig. 10 have strong system size dependence, we 
may conclude that each point in Fig. [TO] might have a (small) 
systematic error. 
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